Thermodynamics of spinless fermion on triangular plaquettes or tetrahedral chain 
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The spinless fermion on triangular plaquettes or tetrahedral chain is considered in this work with 
hooping term, repulsive Coulomb interactions between nearest-neighbor sites and chemical potential. 
This is an exactly solvable model, through decoration transformation and transfer matrix approach, 
when the hopping term acting with nodal sites is ignored, but all other parameters considered are 
quite general. Five different states are shown in the phase diagrams, these states are distinguished 
by particle density, empty, 1, 2, 3 and 4 particles per unit cell on the triangular plaquette and 
one nodal site, or a tetrahedron unit cell. The boundary of these regions are frustrated state 
curves. The phase with particle density 2 is limited by a frustrated state curve with non-trivial 
residual entropy, whereas the other boundaries are limited also by frustrated state curves, but with 
trivial residual entropy. The existence of this non-trivial phenomenon has also been evidenced by 
temperature dependencies of the specific heat. The behavior of the specific heat at low-temperature 
is influenced by the interaction between thermal and quantum fluctuations. Finally some other 
additional quantities were also considered, such as the internal energy and compressibility. 
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I. INTRODUCTION 

The Hubbard model was introduced to describe the 
correlation of electrons in the narrow band of transition 
metals. By extending the model and changing the pa- 
rameters, it is possible that the Hubbard model could 
describe many properties of strongly correlated electron 
system. Montenegro-Filho and Coutinho-Filho [l| con- 
sidered the doped AB2 Hubbard chain, both in the weak 
coupling and the infinite-t/ limit (atomic limit). They 
studied a quite interesting phase as a function of hole 
doping away from half-filled band. In fact, in several 
works, many properties of the system are presented as a 
function of the chemical potential ^. Although from a 
theoretical point of view the chemical potential and par- 
ticle density are conjugated variables, experimentally, it 
is possible to manipulate in a controlled way, by vary- 
ing the entering and outgoing particles to the system. 
Mancini has presented the exact solution of extended 
one-dimensional Hubbard model in the atomic limit, ob- 
taining the chemical potential plateaus of the particle 
density, of the on-site potential at zero temperature and 
studied the thermodynamic charge susceptibility, and the 
compressibility at finite temperature, as well as other 
physical quantities. Vidal et al. Q also discussed two in- 
teracting particles evolving in a diamond chain structure 
embedded in a magnetic field. As the particles interact, 
it leads to the strong localization induced by the mag- 
netic field for the particular value of a flux. Analogous 
model was also studied by Rossler and Mainemer 01 . The 
Hubbard model in other quasi one-dimensional triangu- 
lar structure also was studied by Wang |5|. The latter 
indicated that for small hopping term, the system ex- 
hibits short range antiferromagnetic correlation, whereas, 
when the hopping terms become greater than the criti- 
cal point, there is a transition from an antiferromagnetic 
state to a frustrated one. Moreover, the insulator-metal 
transition takes place at hopping interaction even greater 
than another critical point. Further Gulacsi et al. 0] 
also discussed the diamond Hubbard chain in a magnetic 



field and a wide range of properties such as flat-band fer- 
romagnetism and the correlation induced metallic, half- 
metallic process. 

In the last decade several diamond chain structures 
were discussed. Some of them were motivated by real ma- 
terials such as Cu3(C03)2(OH)2 known as azurite, which 
is an interesting quantum antiferromagnetic model, de- 
scribed by Heisenberg model on generalized diamond 
chain. Honecker et al. Q studied the dynamic and ther- 
modynamic properties for this model. Besides, Pereira 
et al. Q investigated the magnetization plateau of de- 
localized interstitial spins on diamond chain, they also 
detected magnetocaloric effect in kinetically frustrated 
diamond chain. Moreover, the thermodynamics of the 
Ising- Heisenberg model on diamond-like chain was widely 
discussed in Refs. 

EMI, 

and also the thermal entangle- 
ment of this structure was considered . Recently, Lis- 
nii [isj studied a distorted diamond Ising-Hubbard chain 
and that model also exhibits geometrical frustration. The 
spinless versions of the Hubbard model on diamond chain 
was also recently investigated through exact analytical 
solution [l^. Lopes and Dias performed a detailed 
investigation using the exact diagonalization approach. 
The extended Hubbard model on diamond chain struc- 
ture [Tsl also can be solved exactly in quasi-atomic limit, 
and this model exhibits a quite reach phase diagram at 
zero temperature, the thermodynamics also exhibits very 
remarkable properties due to the influence of geometric 
frustration. 

On the other hand, coupled tetrahedral chain struc- 
ture was first introduced by Mambrini et al. [l^ to 
discuss the residual entropy and spin gap in tetrahedral 
coupled Heisenberg chain, later Rojas and Alcaraz (2^ 
also discussed a more rigorous result of the tetrahedral 
Heisenberg chain using numerical approach, the quan- 
tum phase transition and geometric frustration of this 
model, and later Antonosyan et al. [2lj have studied the 
exact solution of the Ising-Heisenberg chain with trian- 
gular XXZ-Heisenberg plaquettes (or tetrahedral chain) 
l2ll. I23 . The Hubbard model version of this structure 
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also was considered by Batista and Shastry [23 , and more 
recently Maksymenko et al. discussed the localized 
states on triangular trap and low temperature properties 
of the antiferromagnetic Heisenberg and repulsive Hub- 
bard models. 

The outline of this paper is as follows. In Sec. 2, we 
introduce the model to be studied. In Sec. 3, we present 
the phase diagram at zero temperature. In Sec. 4 is 
devoted to the decoration transformation and transfer 
matrix approach, in order to obtain its exact solution. 
Whereas in Sec. 5 the thermodynamics properties of the 
model are discussed. Finally, we summarize our conclu- 
sions in Sec. 6. 




FIG. 1: Squematic representation of spinless fermion model 
on triangular plaquettes or tetrahedral chain. The bold lines 
we represent a plaquette with hopping term and Coulomb re- 
pulsion between two sites, while by thin lines we represent 
only the Coulomb repulsion term of nodal particles and par- 
ticles on plaquette site, at i— th site. 



Hamiltonian becomes 



II. TRIANGULAR PLAQUETTES CHAIN 

The one-dimensional lattice of the spinless fermion 
model on triangular plaquettes (2TI . [2^ or also known 
as tetrahedral chain [20| is considered, which is schemat- 
ically illustrated in Fig. [TJ From now on, conveniently 
we will call this structure as the triangular plaquettes 
chain, thus the total Hamiltonian of the model, can be 
expressed by 
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H — y^Ht,i-i_i, 

i=i 



(1) 



where 



Hi,i+i = ^afc^j + al + -aa.i + h.c. 



V 

+-T ("rf.i + ^d,i+i) (na,i + nb,j + n^O , (2) 



and N the number of unit cell; a.a^i{aL^ j) are the Fermi 
annihilation (creation) operators at the ith site, for spin- 
less (or completely polarized) fermion respectively, with 
a = {a, 6, c, c?}, and Wa.i the number operator defined as 

^a,i = i^a,i- The hopping parameter t corresponds to 
the exchange constant of an interstitial spinless particle, 
/i is the chemical potential and the Coulomb repulsion 
term between fermions on neighboring sites are repre- 
sented for V and Vi. 

The spinless fermion model on tetrahedron chain, has 
an important particle- hole symmetry [3, [2^. For the 
purpose of this discussion, the following canonical trans- 
formation j — > a^ i and a^, ^ — > a^ ^ was considered. 
For the occupation number operator the transformation 



— M ^li^fc.i + ^Li^c.'i + a]! jaa,i + /i.i 



+ 2Vi - V) [na,i + rifc^, + n^,,) 

+Vi (na.iHfc.j + n(,,inc^j 4- n^jria^j) 
V 

+ — (rid,,, + nd,i+i) {na.t + nf,,^ -f yIca) 



+ {W - Wi - V) . 



(3) 



Alternatively, using the Wigner- Jordan transforma- 
tion it is possible to map the spinless fermion tetrahe- 
dron chain onto an Ising-Heisenberg triangular plaque- 
ttes chain, the z-component of the spin- 1/2 operator is 
related to the number operator by n = 1/2 — cr^, while 
the creation (annihilation) operators are related through 

a] = 2^ {j{u<, ^l) and a, = 2^ (n,<, -r^) a+. 

In order to solve this model the first step is to obtain 
all the eigenvalues and eigenvectors of the Hamiltonian 
Eq. to evaluate this one, we considered the state 

vector for the sites a, b and c connected by hopping term 
t. In each state there are only two possibilities labeled by 
or 1, that means an empty or occupied particle state 
respectively. 

The states acting on Hamiltonian ([2]) at sites a, h and 
c for each elementary cell become 



1-1,1-1-1 



H 



H 





reads n„ 



n„ , , thus the transformed 



Hereafter, the following notations are introduced, R rep- 
resents a cyclic shift operator, n^; = n^j -|- is the 
occupation number operator for nodal particles of the tri- 
angular plaquettes chain. Finally we define the operator 
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fj- = a± 1 — tR — tR?' , where a+ 
and a- = -2^l + Vi + [V - Hrf. 

After diagonalizing the Hamilton] 
the following eigenvectors 



energies were found, which are described as follows 



V2) 

Va) 

ve) 

Vj) 



I 000), 

-^{1 + R + E?)\ 100), 
v3 
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-1-R+2R') I 100), 



^(-l + i?)|100), 
-^{\ + R + R^)\ 110), 
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-2 + R + R^)\ 110), 



^(i2-i2^)|110), 
I 111). 



(10) 

(11) 
(12) 

(13) 

(14) 
(15) 



The eigenvalues were obtained as a function of the 
chemical potential ^, nearest neighbor Coulomb coupling 
y, Vi and the hopping term t, and the corresponding 
eigenvalues are 



El 

E2 — E3 
Ei 

E5 = Eq 
E7 



A* - 



(16) 
(17) 

(18) 



^2t + Vi-2n+(-^ + V]nd, (19) 



-3/x + 3Fi + 



(20) 
rid. (21) 



where is the eigenvalue of operator Tic; • The energies 
E2 = E3 and i?5 = Eq are twofold degenerate. Further- 
more, using the particle-hole symmetry and considering 
Vi = 1^, we have the following constraints fi 
corresponds to the half-filled band system. 
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III. GROUND STATE PHASE DIAGRAM 

In order to study the phase diagram at T = 0, we use 
the diagonalized version of the Hamiltonian presented in 
the previous section. The ground state phase diagram for 
t as a function of fi is displayed in Fig. ^ and for fixed 
values of Vi = y = 1. In any case, only five ground-state 
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by p we denoted the electron density per elementary cell. 
Therefore, the corresponding eigenvectors are given by 



we obtain 
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(27) 
(28) 
(29) 
(30) 
(31) 
(32) 



where we define the operator g — -^(l + J? + iJ^) and 
the states of type | a, &, c)® | d) corresponds to the par- 
ticles at sites a, 6, c and d, respectively in the unit cell. 

In Fig. [51 it can be observed that there are five ground 
states in the phase diagram. These states are limited as 
follow: 



1 50) : limited by ^ =^ -2t, with p = 0, (33) 
1 51) : in the range -2t ^ /x < 1 with p = 1, (34) 

I j^2) I ■ range 1 < p < 2 with p = 2, 

(35) 

|53) : in the range 2 ^ p < 2t+2> with p = 3, (36) 



1 54) : hmited p ^ 2t + 3 with p = 4. 



(37) 



The first state is |50), this corresponds just to a simple 
empty lattice particle (or fully-filled holes) on triangular 
plaquettes chain with total density p = 0. The second 
one is represented by |51), where one particle may stay 
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FIG. 2: Ground-state phase diagram of triangular plaquettes 
chain, as a function p against t, for V = Vi = 1.0 



at a, b OT c of the elementary cell, then the particle den- 
sity for this state is p = 1. When the density takes the 
value p = 2, we have two-fold degenerate states \F1) and 
\F2), the first degenerate state corresponds to the con- 
figuration in which two particles are fluctuating in sites 
a, b and c. In the second degenerate state \F2), we have 
only one particle fluctuating at sites a, b and c, whereas 
on site d of elementary cell there is another particle. In 
the case of state |53) (p = 3), two particles occupy the 
three interstitial sites a, b and c, in this case, we have 
one particle in nodal site d. Finally there is a fully-filled 
particle in triangular plaquettes chain, whose state can 
represented by |S'4) with density p ~ 4:. 



IV. DECORATION TRANSFORMATION 

In order to study the thermodynamic properties of the 
spinless fermion model on triangular plaquettes chain, it 
will be used the decoration transformation j26l - [29| as de- 
scribed below. It is not necessary to map the spinless 
fermion model onto spin models, like Ising-Heisenberg 
model (ill , [l^ . [30l - l33| . Thus, there is an interest in 
performing the decoration transformation for operators. 
In this sense, it will be applied directly the decoration 
transformation approach [ll|, [l^, [30M33| for the proposed 
model. The main aim to solve the Hamiltonian Eq. ([T]) of 
spinless fermion model on triangular plaquettes chain is 
to map onto an effective spinless fermion model without 
hopping term, whose Hamiltonian is given by 
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Vnd,ind,i+i - - (n<j,i 



(38) 



where V and fl are coefficients to be determined using 
decoration transformation. 

The Boltzmann factor of effective spinless fermion 



model can be expressed as follows 

wind,,, nd,^+l) = e/3(-V'n.,.n,..+ i + #(n,..+„.,.+ i))^ (39) 

where /3 ~ [ksT] ^ with ks as the Boltzmann's constant 
and T as the absolute temperature. On the other hand, 
the Boltzmann factors for the spinless fermion model on 
triangular plaquettes chain is given by the Hamiltonian 
([T|), which reads as 



w (nd,i, n.d,i+i) = tra,b,c (e" 



3Hi 



(40) 



whereas the operators rid.i have two eigenvalues and 
1 , thus, one obtains from Eq. (|40p three different Boltz- 
mann factors, which are written down as 



w{l,l) 



^fiti ^f}(2t+2ti-v) ^ 2e/3(-*+2M-v) ^ 

g/3(2t-yi+3M-2y) _|_ 2g/3(-t-Vi-l-3M-'2y) _|_ 



(42) 



(43) 



For simplicity, we denote by wq = 711(0,0), wi = 
u'(l,0) = ^(0, 1) and W2 = u'(l,l). Applying 
a decoration-interaction transformation, the partition 
function for the present model, must satisfy the following 
condition 



z(F,Fi,/.,t,r) = /Ze//(y,/i,r). 



(44) 



The Boltzmann factor for both systems must be 
equal in order to satisfy the decoration transformation 
(27I. [28j. Thus, there are three unknown algebraic equa- 
tions and three unknown parameters, this algebraic sys- 
tem is solved easily using the decoration transformation 
method, whose solutions are written as 



-M- 



wf 
W2W0 



(45) 



where the /, fl and V are expressed as functions of the 
original parameters of the Hamiltonian by means of wq, 
wi and W2- 

Thereafter, the Boltzmann factor for effective spinless 
fermion model reads as 



Wo 



1, 



Wl 



W2 



(46) 



where x — exp(/3/2/2) and y = cxp(— /3T^). 

We are interested in solving the effective spinless 
fermion without hoping term (also known as in atomic 
limit), to solve this effective model the transfer matrix 
method |34| will be used 



X yx 



(47) 
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It can be noted that the terms ft and V can be ob- 
tained in a similar way as they were obtained for spin 
models. The eigenvalues of transfer matrix are given by 

A± = ^1 + yx^ ± \J{1- yx'^f +4:x'^^ /2. This result is 

exactly the same result to that obtained for the case of 
spinless diamond chain fl6j . 

Using the largest eigenvalues A+ of the transfer matrix 
T, we obtain the partition function per site of the model 
in terms of the effective spinless fermion model in atomic 
limit Z — fZef f , with Zef / = A+ and the partition func- 
tion for effective spinless fermion without hopping term. 
The partition function per elementary cell, is expressed 
in terms of the Boltzmann factors for the spinless fermion 
model on tetrahedron chain, thus the partition function 
becomes 

Z^^(^WQ+W2 + \/(wo - W2)V4w2^ . (48) 

Once known the partition function of spinless fermion 
model on triangular plaquettes chain, it is possible to ob- 
tain the free energy by the relation D, ~ —kT In Z. There- 
fore, we are ready to study several physical amounts such 
as entropy, specific heat, average energy and others. 

V. THERMODYNAMICS 



Before investigating the thermodynamical properties 
of the system, it is worth to observe the chemical poten- 
tial dependence of the total electron density p, as dis- 
played in Fig. [3Ja) for fixed t = 1, V = 1 and Vi = 1, 
and in Fig. for fixed t = 1, y = 1 and V^i = 0.1. At 
very low temperature, the model displays three distinct 
plateaus; for density p = 1, 2 and 3, besides the trivial 
empty particle {p — 0) and fully-filled particle {p = 4) 
states. Here and from now on we will consider all temper- 
ature dependence in unit of fc^ for simplicity For /i = — 2, 
the density jumps from the trivial empty p = to /? = 1, 
where the plateau up to /i = 1. Thereafter, the density 
jumps again to p = 2 and consequently remains roughly 
constant up to /i = 2. The discontinuity for density p ~ 2 
to p = 3 will be present for p = 2. Then the plateau re- 
mains until reaching /i = 5, in which we have a final jump 
of the discontinuity reaching the density p = A. More- 
over, the thermal fluctuation smoothes the curves of par- 
ticle density and for higher temperature just continuous 
crossovers are observed. The curves from different tem- 
peratures also cross roughly at a single point near the 
crossover from the particle density plateau to the satu- 
ration plateau. 

Another important property to be discussed is the sym- 
metry particle-hole described by Eq. (|3]) which will be 
satisfied, therefore, we find the following relation for the 
free energy. 



In this section, we discuss our most interesting re- 
sults, using some relevant thermodynamic quantities that 
enable to reveal the competition between the possible 
ground states. Actually, in several works, many prop- 
erties of the system are presented as a function of the 
chemical potential p. Although from a theoretical point 
of view the chemical potential and particle density are 
conjugated variables, experimentally it is possible to ma- 
nipulate in a controlled way, by varying the entering and 
outgoing particles to the system. Following the same 
procedure to that reference [l^, to obtain explicitly the 
electron density per elementary cell, we take the deriva- 
tive of free energy in relation to the chemical potential of 
spinless fermion on triangular plaquettes chain 



dp 



3(n,) 



(49) 



here we are assuming the exchange invariance of sites a, 
b and c. 

The particle density of type d can be obtained from 
effective spinless fermion model, by the relation 

, , 1 



1 



A+-1 



(50) 



whereas the particle density for sites a, 5, and c, can be 
obtained combining Eqs. p9)) and ((50|) . resulting in 



(51) 



P - (nd> 



n(t, v,Vi) = n{t, V, Vi = v) + (svi - 3V) . (52) 

The chemical potential usually depends on the tempera- 
ture, however at half-filled band electron density p = 2 , 
the chemical potential becomes independent of the tem- 
perature, given simply by /i = ^ . If we look at the phase 
diagram displayed in Fig[3ja), the value p ~ 1.5 corre- 
sponds to the density p = 2. However, the total density of 
electron undergoes true transitions between these phases 
only at zero temperature, where in fact jumps emerge of 
discontinuity in the curves of density in certain critical 
points of the chemical potential (see Fig. Quantum 
fluctuation, rather than thermal fluctuation, are respon- 
sible for such zero temperature phase transitions. 

In Fig. m the temperature dependence of total density 
is shown for different values of the chemical potential p. 
Thus, at r = 0, p leads to discrete values, as shown 
in Fig. m from where two significant situations can be 
explicitly calculated using Eq. firstly, when /i = 1, 

in the limit T — > 0, we obtain p = 2 — while in 
the case p = 2, and T — > 0, the total density leads to 
p = 2 + . These values achieved are slightly different 
from that naive average in an unit cell, which could be 
p = 1.5 or p = 2.5 respectively. This difference is related 
to the frustration and in order to calculate the density p 
we need to calculate the average p in whole lattice instead 
of only an unit cell. Besides, from this plot we may also 
conclude the total particle density {p = 2, for p = 1.5) is 
independent of the temperature only at half-filled band. 
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FIG. 3: (Color Online) Total density p as a function of chemi- 
cal potential /i, for fixed values of the temperature and t = 1.0. 
(a) V = Vi = 1. (b) V = 1 and Vi = 0.1. 
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FIG. 4: (Color Online) Total density p as a function of the 
temperature T, for several values of chemical potential. 



FIG. 5; (Color Online) Temperature dependence of the en- 
tropy for t = 1.0, V — Vi = 1.0 and several values of the 
chemical potential, (a) p = -2.3,-2.1,-2.0,-1.9,-1.7 (b) 
p = 0.7, 0.9, 1.0, 1.1, 1.3. 




FIG. 6: (Color Online) Entropy as a function of the total 
density p for t = 1.0, V = Vi = 1.0 and several values of the 
temperature (a) T = 0.02, 0.06 and 0.1; (b) T = 0.2, 0.5 and 
1.0. 



At this point, we will turn our attention to an inves- 
tigation of the entropy. One can easily determine the 
expression for the entropy according to the general ther- 
modynamic relation 

In Fig. [51 we report the temperature dependence of 
entropy S, for several chemical potential values, for fixed 
t = 1 and V ~ V\ = 1. The behavior of the entropy 
around /x = —2.0 are shown in Fig. El^a), where the 
entropy takes the value In (2) for /i = —2.0, this is the 
so called residual entropy. Moreover, for /i = 5, we have 
the same value of residual entropy. While, in Fig. El (b) 
the residual entropy has a different non-trivial value. To 
obtain explicitly the entropy we return to Eq. (|53p and 
the limit T — ^ 0, the residual entropy results in S = 
In (2) - In (3 - VS) , for /i = 1 and /i = 2. 

At zero temperature, the entropy is discontinuous at 
the transition boundary, signaling the first-order zero- 
temperature phase transition. However, at finite tem- 
peratures, the thermal fluctuations round off the discon- 
tinuity. The total density dependence of entropy are de- 



picted in Fig. IHa). At first, looking more closely in 
the low temperature limit, the degenerate state F\ and 
F2, corresponding to the region of half- filled band parti- 
cle density (p = 2), are broken, as a consequence of the 
entropy increases rapidly. Meanwhile, when the temper- 
ature increases, the peaks are gradually merged into one, 
as shown in Fig. [HJb). 

Another quantity to be plot is the entropy as a func- 
tion of chemical potential, which exhibits a series of peaks 
that correspond to the residual entropy as displayed in 
Fig. [T] Explicitly, the entropy in phase transitions be- 
tween S'O A — > S\ (53 A — !■ 54) takes the value ln(2), 
certainly this factor comes from two different configura- 
tions for each phase, this generates an infinitely degener- 
ated state at zero temperature. However, the transitions 
between 51 A — > F1,F2 {F1,F2 < — > 53) takes a non- 
trivial value of residual entropy 5 = ln(2) — ln(3 — V5), 
this factor again comes from two states, one of them non- 
degenerated and the other two-fold degenerated states, in 
the zero temperature limit. On the other hand, in Fig. 
[7] the behavior of entropy at very low temperatures in 
the neighborhood of the critical values of fi and how the 
residual entropy leads to the value discussed above in the 
limit T = 0. 
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FIG. 7: (Color Online) Entropy as a function of the chemical 
potential jj, with t = 1.0, V = Vi = 1.0, in the low tempera- 
ture limit. 
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FIG. 8: (Color Online) Specific heat C is plotted as a func- 
tion of the temperature T for t = 1.0, V — Vi = 1.0 and 
different values of chemical potential, (a) fi — —2,3, —2.1 — 
2.0, -1.9, -1.7. (b) fj, = 0.7, 0.9, 1.0, 1.1, 1.3. 



A further quantity which is important for a better 
understanding of thermodynamics is the specific heat 
C = ~/3§^, as shown in plots of the specific heat dis- 
played in Fig. [H Thus, for certain values of /i, the spe- 
cific heat has two peaks: one stairdard maximum peak 
and another small peak in the lower temperature region. 
The maximum at low temperature is influenced by the 
phase transition between 5*0, 5*1, Fl, {F2), 53 and 5*4 
states caused by thermal excitations and quairtum fiuc- 
tuations. One can see iir Fig. |S1 a series of peaks that 
appear in the curves at low temperatures, which are in 
agreement with the zero temperature phase trairsition, 
as illustrated in Fig. [51 When the chemical potential 
moves away from the critical point, the peak moves to 
a higher temperature region, roughly keeping its ampli- 
tude. While higher temperatures, the specific heat has a 
different behavior dependiirg on the choice of /i, for ex- 
ample the one in Fig. |8] (a) for fi = —1.9 and /i = —2.1 
at low temperatures the specific heat behavior is similar 
(blue and red solid line), whereas for T > 0.25 the cur- 
vature behavior changes. In addition, for fi — 1.0 and 
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P 

FIG. 9: (Color Online) Specific heat C as function of total 
density p for low values of the T, with t = 1.0 and V = Vi = 
1.0. 



— r=o.oi — r=o.o5 — t=o.\ 

T = 0.2 




12 3 4 

p 

FIG. 10: (Color Online) Total density p dependence of the 
internal energy U, for low values of the temperature T , with 
t = 1.0, V = Vi = 1.0. 

fl ~ 2.0 see Fig. [B^b), the behavior of the specific heat 
becomes remarkably different. 

In Fig. [9] we present the total density dependence of 
the specific heat for three values in the low temperature 
limit, T = 0.02, T = 0.06 and T = 0.1. In this Fig. M the 
corresponding plots exhibit four minimum values which 
are related to the phase transitions at zero temperature. 

In Fig. [TO] we illustrate the behavior of the internal 
energy U = — '^^q^ as a function of the electron density 
p, for different values of the temperature T . In the low 
temperature limit, we can observe several plateaus of in- 
ternal energy that depend on the total particle density. 
The largest interval of almost constant internal energy 
coirfiguratioir occurs between degenerate states | 53) and 
state I 54). Additionally, the smaller jump in the lev- 
els of internal energy configurations occur between states 
|Fl),(|F2))and|53). 

Finally, we conclude our analysis of the thermodynam- 
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FIG. 11: (Color Online) Compressibility k versus total elec- 
tron density p for t = 1.0, V = Vi = 1 and low values of the 
temperature T. 

ics by exploring the compressibility k = ^ for i = 1, 
V = Vi = 1 and different fixed temperatures shown in 
Fig. [TT] The four peaks is a signature of the phase transi- 
tions discussed widely above. The structure of the peaks 
decreases when the temperature increases, in opposite to 
the peaks of entropy (see Figs. Elja) and (b)) increase. 
This means that the system is stiffer close to integer par- 
ticle density, as one can expect, since for those densities 
the system is well organized turning in a rigid structure. 
Due to the particle-hole symmetry, the compressibility is 
invariant under the transformation p — > 4 — p or particle 
hole symmetry. 

VI. CONCLUSIONS 

In the present paper, we studied the spinless fermion 
on triangular plaquettes or tetrahedral chain with hop- 
ping term, repulsive Coulomb interactions particles be- 



tween nearest-neighbor sites and chemical potential (act- 
ing only on triangle plaquette). This model has been ex- 
actly solved using the decoration transformation (26l - [29j 
and the transfer matrix [sJl method. Within the frame- 
work of this exact analytical approach, the ground-state 
phase and the thermodynamics functions of the system 
were discussed in detail. 

The phase diagrams are rather rich, four different 
states and one state two-fold degenerate (|F1) and \F2)) 
are shown, changing from empty lattice of particles up 
to the saturated state of particles on the triangular pla- 
quette or tetrahedral chain. We also observe that, there 
are frustrations between all phase boundaries. The region 
with particle density p = 2 is two-fold degenerated, so the 
boundary of this region becomes a frustrated state curve, 
with non-trivial residual entropy S = ln(2) — ln(3 — \/5), 
whereas the other boundaries also become frustrated 
curves, but with trivial residual entropy S = ln(2), this 
factors comes when two phases becomes identical, pro- 
ducing an infinitely degenerated state, which leads to a 
frustrated state. 

The existence of this non-trivial phenomenon has also 
been evidenced by temperature dependencies of the spe- 
cific heat, as well remarkable characteristic of the 
residual entropy remain constant until relatively high 
temperatures T « 0.3, for chemical potential fi = —2 and 
5. We also discuss the behavior of the particle density as 
a function of the chemical potential and temperature. 

Finally, some other additional quantities were also con- 
sidered, such as the internal energy and compressibility. 
Also we observe that, the behavior of the specific heat at 
low temperature is influenced by the interaction between 
thermal and quantum fluctuations. 
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